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We discuss melting of alloys along grain boundaries as a free boundary problem for two moving 
solid-liquid interfaces. One of them is the melting front and the other is the solidification front. The 
presence of the triple junction plays an important role in controlling the velocity of this process. 
The interfaces strongly interact via the diffusion field in the thin liquid layer between them. In the 
liquid film migration (LFM) mechanism the system chooses a more efficient kinetic path, which is 
controlled by diffusion in the liquid film on relatively short distances. However, only weak coherency 
strain energy is the effective driving force for LFM in the case of melting of one-phase alloys. The 
process with only one melting front would be controlled by the very slow diffusion in the mother 
solid phase on relatively large distances. 



I. INTRODUCTION 

In our previous paper [l| we discuss partial melting 
of the binary alloys which proceeds via the liquid film 
migration (LFM) mechanism. The liquid phase separates 
two nearly parabolic fronts (see Fig. 1). If only the 
melting front existed, the process would be controlled 
by the very slow diffusion in the mother solid phase. In 
the LFM mechanism the system chooses a more efficient 
kinetic path which is controlled by much faster diffusion 
in the liquid film. However, in this case the relatively 
weak coherency strain energy is involved as an effective 
driving force for this process. 

The early observations of liquid film migration (LFM) 
were made during sintering in the presence of liquid phase 
or during partial melting of alloys Q (see [i| for a re- 
view). Nowadays LFM is a well established phenomenon 
of great practical importance. In LFM one crystal is 
melting and another one is solidifying. Both solid-liquid 
interfaces move together with the same velocity. The mi- 
gration velocity is much smaller than the characteristic 
velocity of atomic kinetics at the interfaces. Therefore, 
both solids should at the interfaces be locally in ther- 
modynamic equilibrium with the liquid phase. On the 
other hand, these local equilibrium states should be dif- 
ferent for the two interfaces to provide the driving force 
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FIG. 1: Schematic representation of two moving nearly 
parabolic fronts discussed in Ref. SI and S2 are the 

melting and growing solids, and L is the liquid film. 



for the process. It is by now well accepted (see, for 
example, [J, Q) that the difference of the equilibrium 
states at the melting and solidification fronts is due to 
the coherency strain energy, important only at the melt- 
ing front because of the sharp concentration profile ahead 
the moving melting front (diffusion in the solid phase is 
very slow and the corresponding diffusion length is very 
small) . The solute atoms diffuse ahead of the moving film 
and the coherency strain energy in such frontal diffusion 
zone arises from the solute misfit. Thus, the equilibrium 
liquid composition at the melting front, which depends 
on the coherency strain energy and on the curvature of 
the front, differs from the liquid composition at the un- 
stressed and curved solidification front. This leads to the 
necessary gradient of the concentration across the liquid 
film and the process is controlled by the diffusion in the 
film. 

Thus, a theoretical description of LFM requires the so- 
lution of a free boundary problem for two combined mov- 
ing solid-liquid interfaces with a liquid film in between. In 
Ref. Q this problem was considered for simplified bound- 
ary conditions: the temperature and the chemical com- 
position along each interface were kept constant. This 
means that any capillary, kinetic and crystallographic ef- 
fects at the interfaces were neglected. It was found that 
under these simplified boundary conditions two co-focal 
parabolic fronts can move together with the same ve- 
locity. The situation is rather similar to a steady-state 
motion of one parabolic solidification front into a super- 
cooled melt or one parabolic melting front into a super- 
heated solid. In this approximation the Peclet numbers 
were found, but the steady-state velocity remained un- 
determined at that stage. Thus, the problem of velocity 
selection arises. 

Solvability theory has been very successful in predict- 
ing certain properties of pattern selecting in dendritic 
growth and a number of related phenomena (see, for ex- 
ample, 0, H, [![)■ We note that capillarity is a singular 
perturbation and the anisotropy of the surface energy 
is a prerequisite for the existence of the solution in this 
theory. In [l[ we extended the selection theory for the 
process of liquid film migration where the strong diffu- 
sion interaction between melting and solidification fronts 
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FIG. 2: Schematic representation of the phase diagram a) 
and and configuration of different interfaces near the triple 
junction b). This structure moves steadily with the velocity 
V in the vertical y direction. The triple junction is located at 
the origin. Above the origin two solid grains SI and 5*2 are 
separated by the grain boundary. Below the origin two grains 
are separated by the thin liquid layer L. 



plays a crucial role. 

The nucleation of the melt often takes place at the 
grain boundaries of the mother solid phase. Then the 
thin liquid layer extends along the grain boundary. In 
the present paper we discuss the liquid film migration 
during the partial melting along the grain boundary (see 
Fig. 2). The presence of the triple junction in this geom- 
etry drastically changes the structure of the theory. It 
produces a very strong perturbation of the solid-liquid 
interface and the anisotropy of the surface tension does 
not play an important role in such processes [l(J, EL E2] ■ 



II. FORMULATION OF THE PROBLEM 

We discuss the two-dimensional problem of the steady- 
state motion of a thin liquid film during the process 
of isothermal melting of a binary alloy along the grain 
boundary. The schematic phase diagram of the binary 
alloy and the geometry of the process are presented in 
Fig. 2. We note that the discussed geometry corresponds 
to the spontaneous breaking of the reflection symmetry 
x — > —x. The equivalent structure can be obtained by the 
reflection x — > —x of the structure presented in Fig. 2b. 



We assume that the diffusion in the solid phases is very 
slow and the concentration c in the liquid film obeys the 
Laplace equation. We introduce the normalized concen- 
tration C = (c — cl)/(cl — cs) with cl and cs being 
the liquidus and solidus concentrations of the equilibrium 
phase diagram at a given temperature, To. Then the 
equilibrium concentration and the mass balance condi- 
tions at the solidification front, which separates the grain 
S2 and the liquid phase, read 



Co 



-d K, 



V n = -DdC/dn. 



(1) 



At the melting front, which separates the grain £1 
and the liquid phase, the equilibrium concentration is 
changed by the presence of the elastic coherency strain 
energy Q and also the diffusional flux changes because 
in the solid ahead of the melting front the concentration 
is Co which is different from cs- 

Ci = -bA 2 + d Ki, V n (l - A) = -DdC/dn. (2) 

Here V n is the normal velocity; D is the diffusion coeffi- 
cient in the liquid film; K is the curvature assumed to be 
positive for the interfaces in Fig.l; A = (co— cs)/(cl— cs) 
is the dimensionless driving force; b = Yfl(da / dc) 2 / a 2 f'l 
is the dimensionless constant which describes the co- 
herency strain energy Vt is the atomic volume, Y is 
the bulk elastic modulus, a is the atomic constant, /l(c) 
is the free energy of the liquid phase per atom, f'l is the 
second derivative of /l(c) at c = cj,; do is the chemical 
capillary length which is assumed to be isotropic in the 
present problem, do = 7^V f'l( c L ~ cs) 2 where 7 is the 
surface energy of the solid-liquid interface. 

At an arbitrary rotation of the grain boundary (Fig. 2) 
the angle between fronts at the triple junction 2ip remains 
unchanged, cos ip = 7^/27, where 7;, is the surface tension 
of the grain boundary. This additional thermodynami- 
cal condition eventually allows to find a unique solution 
for the whole structure and the velocity of steady-state 
propagation of this structure. 

In the general case, the formulated problem is quite 
complicated because two fronts interact non-locally 
through the diffusion field inside of the liquid phase. At 
small driving forces A and small angles 2ip between fronts 
this problem can be treated in the so-called "lubrication" 
approximation which is also called the boundary layer 
model (BLM) in the context of solidification. This al- 
lows to reduce the original nonlocal problem to the sys- 
tem of ordinary differential equations for the shapes of 
the fronts. 



III. 



LUBRICATION APPROXIMATION 



The crucial ingredient of the lubrication approxima- 
tion is that the variation of the concentration field in the 
thin liquid layer in the direction normal to the interfaces 
is much larger than in the tangential direction. In this 
approximation the two fronts are close to each other and 
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they can be parametrized by two functions which are in- 
variant to the transformation of the coordinate system: 
K(9) and 8(9). K is the curvature of one of the inter- 
faces and 8 is the distance between the interfaces in the 
normal direction; 9 is the angle between the normal to 
one of interfaces n and the direction of the steady-state 
growth (see Fig. 2) Alternatively, one can use functions 
9(s) and 8(s) where s is the arc-length along the inter- 
face and K = dO/ds. The lubrication approximation is 
valid if 

K8 <C 1 

which is the small parameter of the theory. We will see 
later that this condition is fulfilled if A and ip are small. 

Let us derive the equations for the solidification and 
melting fronts using basic Eqs. (JTJ [3J) and the crucial in- 
gredient of the lubrication approximation that the con- 
centration field varies linearly in the normal direction. 
Then in the main approximation dC/dn w (C\ — C^)/^ 
and from the mass balance condition, Eq. ([T]) wc find 



V cos e = D 



bA 2 - 2d K 



(3) 



where V is the steady-state velocity and we have already 
used the fact that the interfaces are close to each other: 
K\ f» Ki = K. The derivation of the second equation is 
slightly more delicate, because if one takes the difference 
between two mass balance conditions all remaining terms 
are small and one should take into account the small dif- 
ferences between normal velocities and normal gradients 
at two interfaces: 



V nl -V n2 +D 



dC_ 
dn 



dC_ 
dn 



= VA cos ( 



(4) 



Evaluating the term in the right-hand-side which is al- 
ready proportional to the small parameter A wc have 
ignored the difference between normal velocities at two 
interfaces. In the left-hand-side, the small difference in 
normal velocities is due to the differences in the normal 
directions at the two interfaces: 

(Vni - V n2 ) w Vsin 688 /ds = VKsm.6d5/d6. (5) 

The small difference in normal gradients at the two in- 
terfaces is due to the curvature of the interface. The 
simple way to see this effect is to look for the solu- 
tion of the Laplace equation for the concentration field 
C = A+B In r in the local polar coordinate system rather 
than in the Cartesian coordinates as it would be conve- 
nient for flat interfaces. One of the interfaces is located 
at r = 1/K(6) and the other is at r = (l + K8)/K. Then 
using the fact that the gradient scales as 1/r and K 8 <C 1 
we obtain: 
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'd<T 
. dn , 



dC 
dn 



an 



VK8 cos 9. 
(6) 



Inserting Eqs. (|5l[6|) into Eq. ([4| we obtain: 
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K^[8 S m9] 



A cos 9. 



(7) 



Eqs. ([31) and ([7]) are the desired equations in the lead- 
ing order of the lubrication approximation. As seen from 
Eq.© K8 ~ A < 1 which justifies the approximation 
used in the limit of small driving forces. 

At the triple junction, where (5 = and 38/ ds = 2(f, 
Eq. reads: 



tan#n = 



2(p' 



(8) 



where 8q is the value of 9 at the triple junction. The grain 
boundary is rotated by approximately the same angle 9q 
at the triple junction (see Fig. 2). One can eliminate the 
curvature K from Eqs. ([3J[7]) and obtain the closed first 
order differential equation for the function 8(9): 



d_ 



cos ( 



[H(9)sin9] = T —^ 



(0) cos6>' 



where we have introduced the rescaled quantities: 



H = 



SbA 
2d? 



2Vd 
Db 2 A 3 ' 



(9) 



(10) 



We should find the smooth solution of this equation 
which starts from H = at 9 = 9q and diverges as 
H — > l/(^cos0) when 8 — > ir/2. It turns out that this 
is an eigenvalue problem and such a solution exists only 
for specific values of the eigenvalue parameter v which 
depends only on 9$. Using Eqs. ([51 fTU|) we can present 
the final result for the steady-state growth velocity V in 
the form: 



V 



Db 2 A' 
2d 



v(A/2<p), 



(11) 



where the scaling function u(A/2(f), obtained by the nu- 
merical solution of Eq. ©, is presented in Fig. 3 and 
has the following asymptotic behavior: 



and 



v(A/2<p) w 0.25, A/V<1 (12) 



v(A/2<p) « 2.12(A/2^) 3 , A/V>1 (13) 



Eqs. (fTT|) . (fT2|) and ([13]) are the main results of this paper 
providing the expression for the steady-state velocity V 
of the melting process along the grain boundary in terms 
of the driving force A and material parameters. 

The analysis of this section is close in spirit to the 
analysis given in Ref. [l2| where we discussed the melt- 
ing process in eutectic and peritectic systems. However, 
apart from very different physical ingredients involved, 
the lubrication approximations used arc quite different 
in these two cases. In Ref. [l2| the assumption of small 
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FIG. 3: The scaling function v vs. tan#o = A/2ip. 



opening angles tp immediately leads to small variations 
of the angle 9 in the vicinity of the value 7t/2. In the 
present analysis, the additional small parameter A leads 
to wide variations of the angle 9 depending on the ratio 
A/ ip, Eq. ijSJ). Thus, the lubrication approximation de- 
veloped here is more general than the approximation in 
Ref. It involves calculations of the whole scaling 

function v(9q), while in [l2j only a single eigenvalue pa- 
rameter was calculated. Of course, in the limit 9q — > 7r/2 
the present problem can be mapped onto the problem of 
Ref. GJ. 



IV. KINETICS OF THE GRAIN BOUNDARY 

The kinetics of the grain boundary plays an important 
role allowing the necessary rotation of the structure in 
the vicinity of the triple junction ( Fig. 2). The angle 
9 along the grain boundary should change from its value 
9q at the triple junction to tt/2 far ahead of the junc- 
tion. Moreover, the concentration at the grain boundary 
should change from its value c$ at the junction to the 
value Co far ahead. The distribution of the concentration 
c in the grain boundary is given by the solution of the 
Chan equation, [l3| : 



D b 6 b d 2 c/ds 2 - (c - c )V cos 9 = 0, 



(14) 



where D b is the grain boundary diffusion coefficient, 5 b 
is the thickness of the boundary, and s is arc-length; 
s = at the triple junction and s is negative along the 
grain boundary. It turns out that for small driving forces 
A the characteristic length of the concentration decay is 
much shorter than the characteristic length for the shape 
changes. This allows to use the value 9q instead of the 
current values 9(s). Thus, we obtain: 



c(s) - c « (c s - c ) exp(-fc|s|), 



y/V cosOo/ D b 6 b . 



(15) 



The evolution of the grain boundary shape is controlled 
by the recrystallization kinetics [l4[ 



V cos 9 = V b [-d b K b + bA 2 (s) 



(16) 



where V is the steady-sate velocity which we have al- 
ready found from the analysis of the melting process; 
V b is the characteristic velocity scale proportional the 
grain boundary mobility, d b is the capillary length pro- 
portional to the surface tension of the boundary and 
K b = 89 /ds is the curvature. The last term in Eq. 
(Ti"7j]) is due to the coherency strain energy effects and 
A(s) = [c(s) — Co)]/(c£ — cs), where c(s) is given by the 
distribution Eq. (|15|) . Eq. (Til)]) can be integrated once 
leading to: 



(s) -6 = 



Vx(s) bA 2 



V b d b 2kd b 



[exp(-2fc|s|)-l], (17) 



where x(s) = ds cos 9 < (see Fig. 2b). Indeed, the 
characteristic length of the concentration decay scales as 
k^ 1 ~ 1/W and it is much shorter than the character- 
istic length of the angle decay 9(s) which scales as 1/V 
( in the limit of small driving forces A, the steady-state 
velocity V is also small according to Eq. (TTTj) '). We note 
also that the last term in Eq. (jTTj) , which is proportional 
to {D b A/ D) 1 / 2 (kd b ~ (D/D b y/ 2 bA 3 / 2 ), is small in the 
limit of small driving forces A. Thus, far ahead of the 
triple junction the position of the grain boundary relative 
to the junction is given by (see Fig. 2b): 



-(tt/2 - o )V b d b /V. 



Moreover, the diffusional flux along the grain boundary 
has a nonzero value at the triple junction. In principal, 
this flux should be taken into account in the description 
of the diffusional field in the liquid phase. This effect has 
been neglected in the previous sections. More careful 
analysis shows that corrections to Eq. ([3]) and (|7|) due to 
this effect are small in the limit of small A < 1. Thus, 
the grain boundary kinetics plays an important role in 
necessary adjustment of this interface but has no influ- 
ence on the meltin g ki netics as we have already discussed 
previously in Ref. [l2| in the context of the contact melt- 
ing in eutectic and peritectic systems. 



V. DISCUSSION 

As we have already mentioned, the process with only 
one melting front ( see Fig. 4a) would be controlled by 
the very slow diffusion in the mother solid phase. In 
this case the growth velocity scales as in the classical 
dendritic growth. Assuming that the diffusion coefficient 
in the solid phase D s <C D, we can present this scaling 
in the form 



V 



D 2 A* 
Dd Q 



v 7/4 



(18) 
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a) 



b) 



FIG. 4: Schematic representation of the symmetric melting 
structures, a)- with only one melting front during the process 
in bulk of the grain; b)- with only one melting front during the 
process along the grain boundary; c)- combined motion of the 
melting front along the grain boundary and the solidification 
front S3/L via the LFM mechanism. 



This scaling involves the additional small parameter a 
- the anisotropy of surface tension. For melting along 
the grain boundary. Fig. 4b, one should ommit the 
anisotropy factor [111 ] : 



V 



DjA* 

Ddn 



(19) 



With the help of the LFM mechanism the system 
chooses a more efficient kinetic path to relax to the equi- 
librium state. For the case of two nearly parabolic fronts 
without grain boundary (see Fig.l), which was discussed 
in [l[, the migration velocity scales as 



V 



d 



Db 2 A 2 



J/\ A»a 1 / 2 . 



(20) 



In the present paper we discuss melting along the grain 
boundary via the LFM mechanism (Fig. 2b). The nu- 



clcation of the melt takes place at the grain boundaries 
of the mother solid phase. Then the thin liquid layer ex- 
tends along the grain boundary. For small driving forces 
A and for not too small angles (p the steady-state velocity 
of this process scales, according to Eq. (fTTj) . as 



V 



Db 2 A 3 
d 



(21) 



If the nucleation of the melt takes place at the triple 
point where three grains meet together, the process could 
proceed as shown in Fig. 4c. This structure has not been 
investigated theoretically so far. However, we can guess 
that the velocity would scale as in Eq. ([2Tj) . As we 
have already noted, the presence of the triple junction 
in this structure produces a very strong perturbation of 
the solid-liquid interface and the anisotropy of the surface 
tension does not play an important role in such processes. 
Thus, one can use the scaling given by the first relation in 
Eq. (|20|) . which describes the structure in Fig.l, formally 
setting a ~ 1 [ll[ . We note that the velocities of different 
processes increase from Eq. (fT5|) to Eq. (j2"T]) . 



Let us estimate the velocity V from Eq. (f21j) and the 
thickness of the film <5 ~ d n /bA from Eci. tfTD"]) using char- 
acteristic values of the parameters: D ~ 10~ 5 cm 2 /s, 
do ~ 10~ 7 cm, b ~ 0.05 and A - 0.05. This leads to 
V 



10 cm/s and <5 ~ 10 



io- 



In conclusion, we have developed and analyzed a rel- 
atively simple model for the melting kinetics along the 
grain boundary in alloys. The process proceeds via LFM 
mechanism which is controlled by relatively fast diffusion 
in the liquid film. The presence of the triple junction 
plays an important role in controlling the velocity of this 
process. We solved this problem in lubrication approx- 
imation which allows to reduce the originally nonlocal 
problem to a local and analytically tractable problem. 
We derived the expression for the velocity of the melting 
process along the grain boundary in terms of the driving 
force and material parameters. 
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